Graillot et al. SpringerPlus 2014, 3:68 
http://www.springerplus.eom/content/3/1/68 



0 SpringerPlus 

a SpringerOpen Journal 



RESEARCH Open Access 



Coupling groundwater modeling and biological 
indicators for identifying river/aquifer exchanges 

Didier Graillot 1 , Frederic Paran 1 *, Gudrun Bornette 2 , Pierre Marmonier 2 , Christophe Piscart 2,3 and Laurent Cadilhac 4 



Abstract 

Future climate changes and the resulting modifications in anthropogenic activities will alter the interactions 
between rivers and groundwater. The quantification of these hydraulic interactions is absolutely necessary for 
achieving sustainable water use and requires accurate analytical methodologies. This report proposes an 
interdisciplinary approach to the quantitative and qualitative characterization of hydraulic interactions between 
rivers and shallow aquifers, wherein it outlines the advantages of coupling groundwater modeling with biological 
markers. As a first step, we built independent diagnostic maps of hydrological exchanges at the sector scale on the 
basis of hydrogeological modeling and biological indicators. In a second step, these maps were compared to 
provide a quantitative and qualitative understanding of exchanges between groundwater and surface water. This 
comparison significantly improved the calibration of groundwater models through a better assessment of boundary 
zones. Our approach enabled us to identify the conditions under which it could be possible to use biological 
indicators instead of a large set of piezometric measures. The integration of such combined tools in a future 
decision support system will assist governmental authorities in proposing appropriate long-term water policies for 
the preservation of groundwater resources, such as for supplying potable water and/or mitigating pollution risks. 
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Introduction 

The importance of surface water/groundwater exchanges 
that are related to i) river discharges, ii) the availability 
of potable water, iii) river thermal regimes, iv) river 
productivity regarding both nutrient inputs and intersti- 
tial fauna biomass and v) the functioning and dynamics 
of alluvial wetlands has been frequently outlined (Bravard 
et al. 1997; Hayashi and Resenberry 2002; Lucassen et al. 
2005; Jansson et al. 2007). The assessment of groundwater 
fluxes between aquifers, rivers and their floodplains thus 
represents a relevant issue for the long-term management 
of water resources. A precise delineation of the sections of 
aquifers that are mostly supplied by rivers and of river 
sections and floodplain areas that are mostly fed by 
groundwater upwellings (Bornette and Arens 2002) may 
constitute a valuable tool for watershed managers to iden- 
tify strategic water reserves (Griebler et al. 2010). Al- 
though some methods exist for this type of quantification 
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(GICC Research Group 2005), they have two main draw- 
backs. First, building a groundwater model based on a very 
large-scale grid is time consuming. Second, the amount of 
synchronous data that are necessary to feed and fit the 
model is prohibitive. In the framework of global change, 
the preservation of water and the functions that are asso- 
ciated with groundwater/surface water exchanges requires 
new methodologies or combinations of existing methods. 
Such methods should lead to more sustainable water 
usage and, consequently, for the preservation of ground- 
water in good quality and the associated biodiversity 
throughout large rivers. 

The objective of this study was to characterize the hy- 
draulic transfers between a large river and the lateral 
shallow aquifers. Traditionally, in studies of ground- 
water/surface water interactions on the sector scale, hy- 
draulic models have been ubiquitously used to locate 
areas of interactions with aquifers (e.g., the Murray 
River, Lamontagne et al. 2002; the Snake River, Hortness 
and Vidmar 2004; the Jordan River, Holtzman et al. 2005 
and the Loosahatchie River, Urbano et al. 2006). 
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More recently, biological and ecological indicators 
have been used, e.g., microorganisms (Steube et al. 2009; 
Stein et al. 2010) and invertebrates (Schmidt et al. 2007; 
Korbel and Hose 2011). The findings of these studies 
suggest that faunal assemblages reflect hydrological 
exchanges. 

The originality of our study lies in an interdisciplinary 
approach that combines scientific knowledge regarding 
the hydraulic control of groundwater transfers with two 
biological indicators: aquatic plants and groundwater 
fauna, for which the ecological requirements with re- 
spect to groundwater connectivity are known. This ap- 
proach was applied to a large floodplain section to test 
its efficiency in providing the integrative mapping of sur- 
face water/groundwater interactions one a large scale 
(10 to 100-km). 

Materials and methods 

The originality of the proposed method lies in its com- 
bination of physical and biological methods (Bornette 
et al. 2008) which will be successively presented and 
then integrated. 

Site location and description 

The Rhone River runs from Switzerland in the Alps to 
the Mediterranean Sea along a 545-km course. The river 
catchment covers approximately 5220 km 2 of surface 
area (Figure 1). Along this river, several types of aquifers 
occur: free and confined sedimentary aquifers, impervi- 
ous geological formations and aquifers in fissured bed- 
rock. Many man-made structures have been built along 
the Rhone River, including approximately 40 dams for 
hydroelectric power production and flood prevention, 
navigational enhancement, and the management of 
water for irrigation (Bravard et al. 1986; Bravard 1987). 

Our study focused on the Bregnier-Cordon river 
stretch (Ain, France). This site was chosen based on 
three main criteria: i) a low level of human activities; ii) 
the existence of hydraulic and biological data (Creuze 
des Chatelliers 1991; Bornette et al. 1998a) and iii) the 
existence of previous studies (Michal 1988; Bornette 
et al. 1994). This river stretch is located 100-km up- 
stream from the city of Lyon from KP116 to KP91 (KP: 
kilometric distance from a reference point, the Rhone- 
Saone confluence, thus decreasing from upstream to 
downstream; Figure 1). The Rhone River flows from the 
north-east to the south-west, and the slope ranges be- 
tween 0.7 and l%o. The average river discharge is ap- 
proximately 423 m 3 /s, varies from 160 m 3 /s to 1900 m 3 / 
s, and is governed by a pluvio-nival type regime. The 
strong slope is associated with plentiful sediment trans- 
port from the tributaries of the Pre-Alps Mountains, and 
a wide range of braided-style flow regimes has developed 
(Creuze des Chatelliers 1991; Bravard et al. 1986). The 



geology of the area is a complex mix of karstic limestone 
mountains and molassic hills that surrounds a large allu- 
vial floodplain. The left bank of the river consists of par- 
tially karstified limestone formations in the eastern part 
(KP110 to KP100) and a large area molassic hills in the 
western part from KP99 to the downstream end of the 
sector. The right bank consists of karstic limestone and 
marls of the Jurassic period, which is known as Bugey 
Mountain and stretches from KP105 to the downstream 
end of the sector. A small karstic Mountain, Mont de 
Cordon, is isolated between two channels of the river be- 
tween KP2 and KP3 (Figure 1). The floodplain around 
the river is a glacial valley that has been refilled with 
quaternary alluviums. The hydraulic conductivity of 
these alluviums decreases from 1 x 10-3 to 1 x 10-5 m/s 
and the thickness close to the Rhone River is approxi- 
mately 10 m (Chauve 2006). 

The river and its floodplain have been considerably 
modified by a hydroelectric power station that was con- 
structed between 1982 and 1984 (Figure 1). A derivation 
dam divides the river channel into two branches from 
KP104 to KP91 (KP103, Champagneux dam). Upstream 
from this dam, the channel has been transformed into a 
12-km long reservoir (KP112 to KP103). Downstream of 
this dam, the former braided channel (KP103 to KP91) 
is now a by-passed channel with a regulated discharge 
that ranges from 80 m /s in the winter to 150 m /s in 
the summer (Bornette et al. 1994). This by-passed chan- 
nel is also fed by the Guiers River (KP99). Most of the 
Rhone River water is diverted by the upstream dam to 
an 8-km derivation canal (KP1 to KP9) that feeds a 
hydroelectric facility (KP6) with a maximum fall of ap- 
proximately 14.3 m (Q = 780 m 3 /s). To collect water that 
exfiltrates along the banks of the derivation canal and 
the reservoir, two drainage canals have been dug at the 
foot of the reservoir and two at the foot of the derivation 
canal to reduce the potential rise in the floodplain 
groundwater level (Claret et al. 1999). 

Hydraulic diagnosis 

Hydraulic interactions between groundwater and the 
river were preliminarily estimated through gradients that 
were calculated from digitalized piezometric maps with 
GIS spatial analysis tools (ArcGIS 9.2, Arclnfo, ESRI, 
http://www.esri.com/). The method was founded on the 
properties of the triangulated irregular network (TIN) 
model, which is very useful for generating surfaces. 
Each triangle contains three values as defined in Arclnfo: 
elevation, slope (or gradient) and aspect (or angle) 
(Bornette et al. 2008; Paran et al. 2008). 

The proposed method permits the characterization of 
exchanges between any aquifer and river along a linear 
interface (bank) that can be represented as a polyline con- 
sisting of several segments. The exchange characterization 
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Figure 1 Location of the Rhone River in France (top) and map of the study site (bottom). Kilometric points (KP) decrease from upstream to 
downstream for the main and by-passed channels (from KP 1 10 to KP 91) but increase in the deviation canal from the dam (KP1) to the main 
channel (KP 9). 



implies: i) the quantification of the rate of flow (Q in m 3 /j) 
and ii) the identification of the direction of hydraulic inter- 
actions, wherein fluxes from the aquifer to the river are 
characterized by positive values (+Q) and exchanges from 
the river to the aquifer by negative values (-Q). Infiltration 
(Q: m /s) through the bank was calculated using Darcy's 
law (Equation 1) considering gradients (i), the cross- 



section of the river (A: m 2 ) and the permeability of the 
aquifer (K: m/s). 

Q(m 3 /s) = K.A.i (1) 

The gradient values (i) at the aquifer/river interface 
were obtained from digitalized piezometric maps (low 
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water level in February of 1990), which were converted 
in the TIN. Cross-section (A) values were estimated by 
subtracting groundwater elevation values from riverbed 
elevation values. Permeability (K) values were obtained 
from the literature (Michal 1988; Rampnoux 1992; 
Hydratec Design Office 2002). 

The exchange directions were given by the sign (+ or -) 
of the sinus of the angle formed by the groundwater flow 
and river flow. This angle (D) was calculated by subtract- 
ing the angle values given by the TIN of the piezometric 
levels and the TIN of the river levels. The infiltration esti- 
mation (Qe, Equation 2) was weighted by the value of the 
sinus of angle D to take into account the influence of the 
angle between the groundwater and superficial flows on 
the exchanges. The groundwater flow intensity was usually 
too weak when Qe was provided in m /s; hence, Qe was 
given in m 3 /d: 

Qe(m 3 /d) = Q. sin(D). 86400 (2) 

Qe was calculated for each segment of the Rhone river 
bank polyline. The final infiltration estimation (Qf) was 
given by Equation 3. Qf is the sum of each Qe calculated 
for a segment of the polyline multiplied by the length 
(L) of this segment: 

Qf(m 3 /d) =Z (lton) (Qe.L) (3) 

This algorithm facilitated a selection of the scale for 
which the exchange flows had to be calculated. It was 
then possible to aggregate values for several 1-km-long 
unitary segments along a predetermined river Section. 

Interstitial invertebrates 

The interstitial fauna was sampled within the sediment 
using the Bou-Rouch technique (Bou and Rouch 1967; 
Bou 1974), which involves the use of a perforated metal 
pipe. This mobile pipe was driven 50 cm deep at each 
sampling point, where 10 L of a mix of interstitial water 
and sediment was pumped prior to filtering through a 
200-um-size mesh net. Samples were preserved in the 
field in a 5% formaldehyde solution, stained with Bengal 
pink, and sorted in the laboratory. A total of 41 sites 
were sampled in the main active channel of the river, the 
drainage canals and the wetlands located in the flood- 
plain. They were sampled from 1987 to 2009 with vary- 
ing periodicity: 16 stations were sampled in the main 
river channel (sampled twice in 1988 and 1989; Creuze 
des Chatelliers 1991) and three were sampled again in 
2008 and 2009 for the present study. A total of 16 sta- 
tions were investigated in the drainage canals [three sta- 
tions were sampled during the summer in 1987 and 
1995 (Claret et al. 1999), six others during the summer 
in 1995, 1996 and 1997 (Marmonier et al. 2000), and ten 
others in 2008 (unpublished data)]. Finally, five stations 



were sampled in abandoned river channels: two stations 
were sampled in April, July and November in 1987 and 
1995, which were years before and after the restoration 
of the wetland, respectively [fine sediment was removed 
by dredging in 1993 (Claret et al. 1999), and three others 
were sampled in 2007 (unpublished data)]. 

To locate groundwater upwelling zones, we used the 
absolute and relative taxonomic richness (S and S%) and 
the abundances (Q) of the stygobite fauna. We used abun- 
dances and richness because the stygobite fauna is more 
abundant and diverse in upwelling zones than in down- 
welling zones (Marmonier and Dole 1986; Dole-Olivier 
and Marmonier 1992). We also considered relative rich- 
ness values because in areas fed by groundwater inputs, 
the ratio between stygobite and epigean species is ex- 
pected to increase in interstitial assemblages and should 
consequently represent a relevant indicator for hydro- 
logical exchanges. 

Aquatic vegetation 

Aquatic vegetation is strongly governed by the ground- 
water/surface water connectivity in aquatic ecosystems 
(Bornette and Amoros 1991; Bornette and Large 1995; 
Bornette et al. 1998b). The distribution of aquatic plant 
species is related to the nutrient levels (Carbiener et al. 
1990, 1995; Onaindia et al. 1996; Britto et al. 2001) and 
seasonal thermal variation in the water (Bornette and 
Puijalon 2009, 2011). In aquatic habitats, there is a 
strong correlation between the intensity of groundwater 
exfiltrations and seasonal temperature variability (Sakura 
1993; Bornette and Large 1995). The various water 
sources that supply aquatic habitats usually strongly 
contrast in terms of nutrient levels (Bornette et al. 
1998a; Bornette and Arens 2002). Knowledge of the eco- 
logical requirements of species that occur in a given 
habitat allows the determination of the origin and turn- 
over rate of the water that supplies the habitat (Amoros 
et al. 2000). For this purpose we had to i) identify the 
ecological range of aquatic species in terms of trophic 
levels in the water and temperature, ii) map the occur- 
rence of aquatic species in the aquatic habitats of the 
floodplain and iii) determine the weight of groundwater 
supplies in each of these aquatic habitats and their po- 
tential origin, as estimated using the abundance and eco- 
logical requirements of the aquatic species. 

Identification of the ecological range of aquatic species 

The ecological range of each species was based on the 
occurrence of a significant relationship between species 
occurrence and relative abundance and i) the thermal re- 
gime and ii) the water trophic levels of wetlands. For this 
purpose, we used floristic releves that were collected 
from 1993 and 1999 in 36 wetlands of the Rhone and 
Ain Rivers. The physico-chemical characteristics of the 
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wetlands were simultaneously measured monthly for at 
least a year in the same set of wetlands (i.e., 12 replicates in 
the Rhone River wetlands and 24 for the Ain wetlands in 
terms of ammonia nitrogen, phosphates and temperature). 
This sampling was performed in the same year as (Rhone 
River) or the same year and year after (Ain River) vegeta- 
tion sampling. At the same time, three to six substrate 
samples were collected in all of the wetlands to determine 
the phosphorous level of the substrate. Ecological profiles 
were produced for all species occurring in at least five sites 
among the 36. 

Cool, well-oxygenated groundwater supplies limit 
plant growth and competition processes but also allow 
evergreen species to grow throughout the year (Bornette 
and Large 1995; Greulich and Bornette 2003). Conse- 
quently, the distribution of aquatic plants is partly ruled 
by water temperature. To quantify the groundwater sup- 
plies to the wetlands, the relative standard deviation of 
temperatures was calculated for each wetland, and the 
wetlands were divided into five equilibrated classes ac- 
cording to this coefficient. In a second step, the average 
abundance of each species for all of the wetlands belong- 
ing to a given class was calculated. The difference be- 
tween the classes was then tested for each species. First, 
the equality of variances between classes was tested 
using a Levene test. Second, either parametric or non- 
parametric tests were performed to determine whether 
the abundance of species significantly differed between 
classes. 

Combined with knowledge about groundwater supplies, 
the nutrient content of the water in the wetlands provides 
information about the origin of groundwater, that is, the 
seepage water is nutrient rich in the Bregnier-Cordon area, 
whereas the karstic lateral aquifers should be less eu- 
trophic (Bornette et al. 1994, 1998a). First, the species af- 
finity for phosphorous (both in the water [PO4 ] and the 
substrate total P content) was assessed. For this purpose, 
monthly [PO4] values were log transformed and averaged 
per site, and the total phosphorous values were standard- 
ized and log transformed. Then, both values were summed 
for each site, and the same method that was used for 
temperature was applied: wetlands were categorized into 
five classes according to their total phosphorous contents, 
and the non-random character of the species distribution 
along the phosphorous gradient was tested. The same pro- 
cedure was applied to the ammonia nitrogen contents of 
water [N-NH 3 ]; however, wetlands were separated into 
four instead of five equilibrated classes. The species affin- 
ity for nitrate nitrogen is not assessed because it has been 
demonstrated to rarely rely to eutrophication. Phosphor- 
ous is frequently the major limiting nutrient for aquatic 
plants (Amoros et al. 2000; Bornette and Puijalon 2009, 
2011). The affinities of aquatic species for phosphorus and 
ammonia nitrogen were then qualitatively combined in a 



trophic profile of four classes by evaluating the relative af- 
finity of the species for eutrophication. The analyses were 
carried out with Jump (JMP* Statistical Discovery Soft- 
ware, http://www.jmp.com/). 

The mapping of the aquatic species occurrence 

The aquatic species contents and relative abundances 
were then roughly determined for each aquatic habitat 
of the Bregnier-Cordon reach. For this purpose, at least 
one floristic releve was carried out in each wetland in an 
area that was considered to be representative of the flor- 
istic contents of each wetland. For rivers, canals and 
riverine wetlands, several releves were performed to take 
into account changes in water quality along the up- 
stream/downstream gradients. For each releve, the rela- 
tive abundance of each species was estimated using the 
Braun-Blanquet (1932) index. Then, the trophic levels 
were estimated by i) multiplying the abundance of each 
species by its ecological profile and ii) summing the 
values obtained for all species of the releve (Amoros 
et al. 2000). We thus obtained the affinity of the wetland 
for each trophic class. We considered the class present- 
ing the highest values to indicate the most probable 
trophic level (Additional file 1). For groundwater inputs, 
each species was classified according to its affinity for 
stenothermic water. The difference between the number 
of species that were intolerant to stenothermic waters 
and those preferentially occurring in groundwater-fed 
ecosystems was determined. Then, the wetlands were 
mapped according to four groundwater influence classes 
(null: index value = 0; low: index value < 4; intermediate: 
index values ranging from 5 to 9; high: index value > 9) 
and four trophic level classes (oligotrophic, mesotrophic, 
eutrophic and hyper-eutrophic). A comparison of the two 
maps then facilitated the identification of several geo- 
graphic areas with contrasting groundwater influences. 

Map crossed analysis 

The three types of describers resulted in the production 
of three separate sets of maps that all located the shal- 
low aquifer recharge areas by the river and sections of 
the river fed by groundwater. These maps cannot be dir- 
ectly compared because of their characteristics; however, 
it was possible to analyze the global patterns that re- 
sulted from the three different describers using overlay 
mapping. 

A this stage of our work, we can only analyze the glo- 
bal patterns because the available data were not col- 
lected simultaneously. However, the crossed analysis is 
possible because data were sampled during comparable 
hydrologic conditions (low water level), during the same 
decade (February 1990 for piezometric values; summers 
1987, 1995, 1996, 1997 and 2008 for interstitial in- 
vertebrates data; summers 1993 to 1999 for aquatic 
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vegetation and chemical data) and after the modifica- 
tions inferred by the hydroelectric power station. Fur- 
thermore, we will look if significant correlations can be 
found between estimated infiltration, stygobite richness 
and abundance and macrophyte index. 

Results 

The diagnoses of the surface water/groundwater ex- 
changes that were made using the three describers were 
first independently considered for the three types of 
metrics, and the subsequent three maps were then 
superimposed. 

Hydraulic modeling 

Gradients (i) were obtained from the piezometric map 
TIN (low water level in February of 1990). The average 
cross-section of the river (A) was approximately 4.2 ± 
1.2 m (min = 1.8, max = 8.1). The following permeability 



values (K) were obtained: Area 1, 1.3 x 10-3 m/s; Area 
2, 5x10-3 m/s; Area 3, KP111 to KP106, K = 5xl0- 
3 m/s, KP106 to KP102.5, K = 3.9x10-4 m/s and 
KP102.5 to KP99, K = 8.5 x 10-3 m/s (Figure 2). The in- 
filtration estimation (Qe) was calculated for each area. 
For each of the zones with the same global exchange be- 
havior, we calculated the Qf and the unit flow (Qu) for 
1-km (m 3 /day/km). 

From a hydrogeological point of view, the GIS spatial 
analysis method was tested on three distinct aquifers 
to calculate groundwater/river exchanges (Figure 2 and 
Table 1): 

- First, the polyline of Area 1 (the La Charriere zone, 
approximately 5600 m long, KP116 to KP112) 
consisted of 62 segments. The results provide 
evidence that the left bank of the Rhone River can 
be divided into two zones (Figure 2): zone 1 (KP116 
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Figure 2 Mapping of the piezometric levels in the alluvial aquifer and the hydraulic exchanges between the river and the 
alluvial aquifer. 
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Table 1 Synthesis of the exchange flow between groundwater and Rhone River in the Bregnier-Cordon sector 


Name 


Stretch / KP 


Flow direction 


Qf 

m 3 /day 


Bank length m 


Qu 

m 3 /j/km 


Area 1 (Al) 


Zone 1 (Z1) 
1 16 to1 14.5 


- 


-668 


1547 


-432 


La Charriere 


114.5 to 112 
Zone 2 (Z2) 
Total Al (11 +Z2) 
_ 


+ 


-52 

393 

-720 


531 
3545 

2079 


-98 
111 

-347 




+ 




393 


3545 


111 


Area 2 (A2) 


Zone 3 (Z3) 
1 1 3 to 1 1 1 




-1890 


2040 


-927 


Bovinel 


Zone 4 (Z4) 
1 1 1 to 1 09 


+ 


6097 


2809 


2170 


Area 3 (A3) 


Zone 5 (Z5) 


_ 


-354 


587 


-603 




1 1 0.5 to 1 1 0 


+ 


5 


76 


64 




Zone 6 (Z6) 




-77 


431 


-178 




1 1 0 to 1 06 


+ 


1723 


3709 


465 


Les Rubattiers-Leschaux 


Zone 7 (Z7) 




-462 


3057 


-151 




1 06 to 1 02.5 


+ 


8 


349 


24 




Zone 8 (Z8) 


- 


-35 


21 


-1693 




102.5 to 99 


+ 


54340 


4190 


12968 


Total A3 (Z5 + Z6+Z7 + Z8) 






-928 


4096 


-226 


+ 






56077 


8323 


6737 


Total A1 + A2 + A3 






-3538 


8214 


-431 


+ 






62566 


14678 


4263 



The zones (Z) and stretches (S) are shown in Figure 2. 



to KP114.5), which is associated with the river 
supplying the aquifer and a Qf of -668 m /day, and 
zone 2 (KP114.5 to KP112), which consists of an 
aquifer drained by the Rhone River and a Qf of 
393 m 3 /day. 

- Second, the polyline of Area 2 (the Bovinel zone, 
approximately 4850 m long, KP113 to KP109) was 
comprised of 74 segments. The results provide 
evidence that the left bank of the Rhone River to 
can divided into two zones (Figure 2): zone 3 
(KP113 to KP111), with the river supplying the 
aquifer and a Qf of -1890 m /day, and zone 4 
(KP111 to KP109), with an aquifer drained by the 
river and a Qf of 6097 m 3 /day. 

- Third, the polyline of Area 3 (the Les Rubattiers- 
Leschaux zone, approximately 12400 m, KP111 to 
KP99) consisted of 312 segments. These results 
provide evidence that the left bank of the Rhone 
River can be divided into four zones (Figure 2): zone 



5 (KP111 to KP110), with the river supplying the 
aquifer and a Qf of -354 m 3 /day; zone 6 (KP110 to 
KP106), with the river draining the aquifer and a 
Qf of 1723 m 3 /day; zone 7 (KP106 to KP102.5), 
with the river supplying the aquifer and a Qf 
of -462 m 3 /day; and, finally, zone 8 (KP102.5 to 
KP99), with the river draining the aquifer and a Qf 
of 54340 m 3 /day. 

Overall, the Rhone River in the Bregnier-Cordon sec- 
tor gains 62566 m 3 /day of water (~4263 m 3 /day/km) 
and loses 3538 m /day of water (-431 m /day/km). A 
general finding which is not explicitly mentioned is that 
the river gained much more water from the aquifer than 
it lost to the aquifer within the stretch investigated. In- 
deed, the groundwater discharge into the river was very 
low (i.e., ~0.72 m 3 /s) in comparison to the Rhone River 
discharge (i.e., ~160 m /s during the season of low 
discharge). 
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Table 2 List of stygobite species sampled in the Bregnier-Cordon sector between 1988 and 2008 


Order 


Species 


KP of occurrence 


Mesogastropoda 


Islamia moquiniana 


3 - 6; 97 - 1 03 


Ostracoda 


Fabaeformiscandona wegelini 


97 - 1 03 


Bathynellacea 


Parabathynella cf stygia. 


96 - 97 


Amphipoda 


Crangonyx subterraneus 


96 - 97 




Nipgargopsis casparyi 


4 -5; 96 - 103; 109- 110 




Niphargus fontanus 


4 - 5; 96 - 98; 109 - 110 




Niphargus kochianus 


4 - 6; 97 




Niphargus rhenorhodanensis 


96 - 103; 109 - 110 




Salentinella juberthiae 


96 - 97; 109 - 110 



Stygobite fauna 

A total of 106 taxa were sampled, of which nine were 
stygobite species (Table 2). The local richness, abun- 
dance and composition of the stygofauna greatly varied 
with the channel type and location (Figure 3). 

In the by-passed channel, two contrasting sections were 
distinguished. From the dam to the confluence with the 
Guiers River and to the beginning of the anastomosed sec- 
tor (KP103 to KP96), seven stations harbored stygobite 
species in their interstitial assemblages (Section a, Figure 3). 
The abundance and richness of stygobite species greatly 
varied, with locally significant values being obtained, espe- 
cially for the varied abundance (15 individuals at KP101 



and 13 individuals at KP99), suggesting punctual inputs of 
groundwater into the by-passed channel. Here, the stygo- 
bite species consisted of shallow interstitial fauna (Fabae- 
formiscandona wegelini, Islamia moquiniana, Niphargopsis 
casparyi and Niphargus rhenorhodanensis). Downstream 
of this section (KP91 to KP96), no stygobite fauna was 
sampled in the by-passed channel, suggesting little or no 
groundwater inputs there (Section b, Figure 3). This result 
was consistent with the geomorphological characteristics 
of this section, which is associated with sediment depos- 
ition and clogging. 

Additionally, contrasting results were obtained in the 
four drainage canals (Figure 3). Upstream of the dam 
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(KP112 to KP109), the abundance and richness were 
high at the upstream end of the right and left side drain- 
age canals, suggesting important groundwater inputs 
(Section c, Figure 3). Here, the stygobite fauna consisted 
of two crustacean amphipod species living in shallow 
interstitial habitats of the main channel (N. casparyi, N. 
rhenorhodanensis) and two species living deeper in the 
sediment (Salentilella juberthiae), or frequently found 
far from the main channel (Niphargus fontanus). This 
last species may indicate inputs of groundwater from the 
nearby karstic system on the left bank. From KP109 to 
KP103, in the area influenced by the reservoir, no stygo- 
bites were sampled in the drainage canals, suggesting a 
weak influence of groundwater in comparison to import- 
ant inputs of surface water that filtrate through the river 
banks (Section d, Figure 3). In contrast, in the left and 
right side drainage canals that are located downstream 
of the dam (KP4 to KP6), stygobite species were sampled 
in the entire length of the left canal and at the upstream 
end of the right one, but with low richness and low 
abundances (Section e, Figure 3). Here, the stygobite 
fauna consisted of two amphipod species living in shal- 
low interstitial habitats of the main channel (N. casparyi, 
N. rhenorhodanensis) and four species living deep in the 
sediment far from the main channel (N. fontanus, 
Niphargus kochianus, Crangonyx subterraneus and Para- 
bathynella cf stygia), suggesting a mix of surface water 
infiltrating through the bank and deep groundwater 
coming from the nearby karstic systems (i.e., Bugey for 
the right-side canal and Mont Cordon for the left one). 

Aquatic vegetation 

The observed aquatic vegetation facilitated the discrim- 
ination of four areas. In the first upstream area, the 
Rhone River flowed as a unique channel along a karstic 
massif. All of the wetlands were oligotrophic or meso- 
trophic, with the groundwater influence ranging from 
low to high. No eutrophic wetlands occurred along the 
right side of the river, suggesting that the wetlands were 
not supplied by seepage water but, instead, by a drained 
hillslope aquifer in this area. Along the left bank of the 
river, some hyper-eutrophic wetlands occurred along the 
karstic massif, suggesting either that the wetlands did 
not drain the hillslope aquifer but, instead, were supplied 
by river seepage, or that the aquifer was eutrophic in this 
area. Downstream of this reach, a zone where the wet- 
lands appeared to be eutrophic occurred. Because the 
river seepage is hyper-eutrophic in this area (Bornette 
et al. 1994), this pattern could indicate either the drain- 
age of a eutrophic aquifer by the wetlands, or that the 
wetlands are supplied by a mixture of both oligotrophic 
groundwater and hyper-eutrophic seepage. This second 
interpretation was reinforced by the occurrence of oligo- 
trophic wetlands along the floodplain margins. A third 



zone was discriminated by the braided zone of the river, 
wherein the sampled vegetation indicated that most wet- 
lands were poorly supplied by groundwater and were 
hyper-eutrophic, indicating that the river was aggrading 
in comparison to its floodplain, and that only a few 
seepage-supplied wetlands occurred. Along the river 
margin, the wetlands were highly supplied by ground- 
water and less eutrophic, suggesting some influence of a 
hillslope aquifer. 

Superimposition of the maps 

The superimposition of the three maps (Figure 4 and 
Table 3) highlighted the convergences and divergences 
among the describers: 

- In most cases, the results of the three methods were 
consistent. In area 1, at La Charriere (left bank, 
KP114.5 to KP112), there was a lack of groundwater 
invertebrates; however, hydraulic modeling and the 
macrophyte index both indicated poor to low (Qu = 
111 m 3 /day/km) inputs of groundwater to the river. 
From KP110 to KP106 (left bank), the macrophyte 
index and the hydraulic method converged and 
indicated low intensity drainage of the aquifer by the 
river (Qu = 465 m 3 /day/km). From KP106 to 
KP102.5 (on the left bank), the three methods 
provided similar results, i.e., a low supply of river 
water to the aquifer (Qu = -151 m 3 /day/km). Finally, 
the results were also consistent for the KP102.5 to 
KP99 section (left bank) and indicated medium to 
high intensity drainage of the aquifer by the river 
(Qu = 12968 m 3 /day/km). For the right bank, all of 
the methods provided similar results for the reach 
located from KP110.75 to KP109, i.e., medium to 
high intensity drainage of the aquifer by the river 
(Qu = 2170 m 3 /day/km). 

- Discrepancies only occurred between the different 
methods in two cases. First, in zone 2, at Bovinel (right 
bank, KP113 to KP111), the macrophyte index did not 
detect a flow from the river to the aquifer, whereas the 
hydraulic method did (Qu = -927 m 3 /day/km). 
Second, in zone 3, at Les Rubattiers-Leschaux (left 
bank, KP111 to KP110), the three methods diverged 
over a small distance: the macrophyte index did not 
highlight groundwater inputs, whereas the hydraulic 
method indicated a weak flow from the river to the 
aquifer (Qu = -603 m /day/km) and the invertebrate 
index indicated an important mixing of surface water 
and groundwater in the drainage canal. 

- Finally, in the zones where no hydraulic diagnosis 
was available (downstream of KP99 and along the 
derivation canal), the macrophyte and groundwater 
fauna indices gave rather similar results. From KP4 
to KP9 (right bank, derivation canal) and around 
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Mont de Cordon (left bank, derivation canal), both 
biological indices indicated groundwater upwelling 
to the drainage canals and to the river. Downstream 
of KP99 (left bank) and KP9 (right bank) to the end 
of the study area, both macrophyte and invertebrate 
indices indicated a lack of groundwater influences. 

Significant correlations were observed between the 
three methods (i.e., estimated infiltration (Qu), stygobite 



richness and abundance and macrophyte index, n= 11 
zones). Thus, we found a significant correlation either 
between the macrophyte index and invertebrate rich- 
ness (R = 0.407, p = 0.034) and invertebrate abundance 
(R 2 = 0.382, p = 0.042) or between the macrophyte index 
and the infiltration estimation (R 2 = 0.503, p = 0.014). 
No significant correlation was found between the hy- 
draulic infiltration estimation and the groundwater 
invertebrates. 




Legend 



A KP 

Hydraulic Characterization 

Rhone River inputs (Qu in m 3 /day/km) 

Medium (- 1284 / -583) 

- Weak (- 157/- 151) 



Aquifer inputs (Qu in m /day/km) 
Weak (111 /465) 

Medium (2170) 

High (12968) 

Figure 4 Comparison map that combines the results from the hydraulic, invertebrate and macrophyte indicators. 
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Table 3 Comparison of the results of the hydraulic and biological methods 

Hydraulic 



Qf Qu Macrophytes Invertebrates 



KP 


Bank 


Type 


m 3 /day 


m 3 /day/km 






115J14.5 


Right 


Rhone River 






Weak, karst 




114.5J13 


Right 


Rhone River 


: 


: 


Nil 




1 1 3_1 1 1 


Right 


Rhone River 


-1541 


-871 


Nil 


- 


111 1 1 0 75 


Right 


Plam infh ipnrp 

L^Ctl II II II I U til IL.C 


-349 


-1 294 


K/loHiiim tn hinh L'arcit 
1 VI CU 1 Ul 1 1 LU 1 1 iy 1 1, i\a\ j L 




110 75 109 

1 IW,/ J 1 W J 


Right 


Pis m i nf 1 1 ipnrp 
LyQi 1 1 ii iiiutri iL_cr 


6097 


21 70 


K/lpHiiim tn hinh kar^t 

IVItTLJILJIII tU 1 1 1 1 1 , i\Q 1 J L 




1 09 1 08 


Right 


Dsm infh ipnrp 
LyQi m iiiiiutriiL_tr 






MpHiiim tn hinh kar^t 

1 V 1 tTLJ 1 LJ 1 1 1 tU 1 IIUL 1^ i\0 1 J L 




1 08 1 05 


Right 


Plain inf i ipnrp 
L/al II II 1 1 1 Utrl IL.C 






Mi tn \rt/pal" karct" 
inii lu vvtrdtx, r\diii 




1 05 1 03 


Right 


Dsm infh ipnrp 

VJO. 1 l 1 IIIIIUtrilL_tr 






\A/pal<- tn mprlii im kar^t 

V v Cul\ LU 1 1 1 tTLJ 1 U 1 1 1, r\u 1 J L 




1 4 


Right 


DerivBtion cbfibI 










4 9 


Right 


Pipriv/atinn rp\ n?i 1 

L^tTI 1 V u I IW 1 1 L.C3 1 lul 

Restitution 






K/lpHii im tn hinh ^ni lifpr nf thp f-ilanrl 

fVltTtJIUIII l\J. 1 IIUL I, QLJ Li J 1 CI U LI 1 LT LJ 1 G 1 1 LJ 

River 


nppn n rni i nrlw^tpr flnw 

LVLTLTL^ LJIWUIILJVvuLLTI i ilj vv 


9 Hnw^trpsm 

J UUVVjLICOI I 1 


Right 


Rpctjti itinn Rhnnp 
i \tr _) l i l u 1 1 <j i i ; im iui itr 

River 






VA/p^k tn nil 

v V L7G r\ LU 1 1 1 1 


Wpal nrni inrlwstpr fln\A/ 


116 115 


Left 


Rhnnp Rivpr 
iMiuiitr mvtri 


-582 


-583 






115J14.5 


Left 


Rhone River 


-86 


-157 


Nil 


- 


114.5J12 


Left 


Rhone River 


393 


111 


Nil 




1 1 2 J 1 1 


Left 


Rhone River 






Nil 




111J10 


Left 


Dam influence 


-354 


-603 


Nil, pollution from karst 


Important exchanges, mixing 


1 1 0_1 06 


Left 


Dam influence 


1723 


465 


Nil to weak 




106J02.5 


Left 


Dam influence 


-462 


-151 


Nil to weak, flow from the Rhone River 
to the aquifer 


Flow from the Rhone River to 
the aquifer 


102.5_99 


Left 


Rhone River 


54340 


12968 


Medium to high, groundwater 


Flow from the aquifer to the 
Rhone River 


99_dowstream 


Left 


Rhone River 






Weak, preferential corridor, groundwater 


preferential corridor, 
groundwater 


103_99/ 1_5 


Mont 
Cordon 


Rhone River, 
Derivation canal 






Karstic reappearance 


Karstic reappearance 



Discussion 

Validation elements for hydraulic exchanges 

Hydraulic studies have been previously conducted to es- 
timate exchanges between the Rhone River and ground- 
water in the same area of Bregnier-Cordon based on 
hydraulic modeling and different piezometric measure- 
ments (Rampnoux 1992; Michal 1988; Kosmelj 1982) 
but without the quantification of fluxes. These studies 
provided hydraulic gradient values for five zones along 
the right and left banks of the river. The directions of 
these exchanges were confirmed by the results of the 
present study. The quantifications of the fluxes that were 
performed in the present analyses highlighted very weak 
values that could be associated with groundwater fluxes. 
This suggests that the Bregnier-Cordon aquifer repre- 
sents a limited reserve of groundwater and that its use 
for irrigation or as a source of potable water must be re- 
stricted. The multidisciplinary strategy employed in this 
study could lead to a new modeling approach that can 



be applied on a broad scale in instances where it would 
be too costly and time consuming to implement 
methods that use gridded datasets. 

Multi-metric comparison 
Convergence of the results 

Despite the absence of significant correlation between 
hydrology and groundwater fauna, we observed a sur- 
prisingly good convergence of the results that were ob- 
tained using hydraulic and biological indicators. This 
was especially observable on the left bank of the river 
from KP106 to KP102.5. In this section, which was influ- 
enced by the dam, the estimated flow value was 460 m 3 / 
day from the river to the aquifer. The value of the 
macrophyte index was very low and indicated a lack of 
water exchange between groundwater and surface water. 
The invertebrate index confirmed the direction of the 
flow and its intensity in comparison to other sections of 
the river. Between the dam and the confluence of the 
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Guiers River at KP99, the flow rate value significantly in- 
creased until reaching 45.367 m 3 /day from the aquifer to 
the river. This result was confirmed by both the macro- 
phyte and invertebrate indices. 

Divergence needing complementary approach 

Interestingly, we did not find this convergence at all of 
the sites. For example, between KP111 and KP110, the 
hydraulic gradient and macrophyte index demonstrated 
that the river fed the aquifer; however, the invertebrate 
indicators highlighted important water exchanges. In this 
case, further investigations through geochemical and iso- 
topic analyses may facilitate the identification of the ori- 
gin of the water. The use of environmental isotopic 
indicators would be a valuable method for validating the 
results of the present study (Fette et al. 2005). For ex- 
ample, the ratios of isotopic hydrogen (5D) or oxygen 
(5 ls O) are usually used for identifying river/aquifer ex- 
changes (Gonfiantini et al. 1995; Walker and Krabbenhoft 
1995). This approach has been used on other sections of 
the Rhone River to identify the origin of surface-water 
and groundwater (Paran et al. 2012). 

Conclusion about multi-metric comparison 

Moreover, a combination of these three methods pro- 
vides complementary information (Bornette et al. 1998a; 
Bornette and Arens 2002; Paran et al. 2005). The hy- 
draulic method precisely quantifies the flow and direc- 
tion of the exchange, whereas the macrophyte and 
groundwater invertebrate indices may reveal the relative 
importance of groundwater inputs in comparison to sur- 
face water flows and can help to distinguish the origin of 
groundwater. 

Method limitations and accuracy 

The accuracy of the results provided by the method dis- 
cussed herein is different for each metric. For hydrogeol- 
ogy, we find the traditional sources of error, such as the 
depth of the aquifer, discontinuities in the alluvial aqui- 
fer and measurement errors related to groundwater level 
and hydraulic conductivity. For hydraulics, the depth of 
the channels has been estimated to be up to 5 m based 
on the quality of the data provided by the National 
Company of Rhone. The possible clogging of the banks 
has been estimated from the results given by the models 
when they were available. 

For invertebrates, the limitations of the method lie in: 
i) the time required for sampling, ii) the identification of 
species, which is also time consuming and requires a 
large panel of specialists; and finally, iii) the lack of 
knowledge about the ecological preferences of many 
groundwater species. For the Bregnier-Cordon sector 
and for most of the northern part of Europe and North 
America, another limit related to the use of groundwater 



fauna is their low diversity in biogeographic zones that 
were frozen during the last glaciations (Ferreira et al. 
2007). Most populations of groundwater species disap- 
peared under quaternary glaciers and recolonized these 
areas during the last 10 ky (Ferreira et al. 2007). This 
low species richness certainly reduces the efficiency of 
the groundwater invertebrate index. Future research 
must test this index in areas located outside of the qua- 
ternary glacial limit. 

Aquatic vegetation surveys are not time consuming in 
comparison to analyses of aquatic invertebrates, as most 
species can be identified in the field. Furthermore, the 
ecological requirements of plants are usually rather nar- 
row because of their sessile character, improving their 
indicative value. However, there are some limitations re- 
lated to i) the paucity of wetlands in the floodplain and 
ii) the existence of poorly informed species, either be- 
cause they are rare and not documented, or because they 
exhibit a broad ecological niche. Finally, it is still neces- 
sary to test whether these ecological profiles are usable 
for other river systems. 

At this stage of our work, we know that our method 
works but remains perfectible. To improve the method 
and its accuracy, our future works will focus on syn- 
chronous data to reduce the risk of errors in the inter- 
pretation of results and to make the statistical analysis 
efficient. Furthermore, we have to better consider the 
scale of interpretation: the spatial resolution of the esti- 
mated infiltration and the macrophyte index which is 
about 1-km versus 100 m for the stygobite richness and 
abundance. 

Conclusions 

The numerical modeling of groundwater/ surface water 
interactions is not always efficient for large-scale river 
basins. In very well-instrumented areas, the exchanges 
between groundwater and surface water can be calcu- 
lated from a spatial analysis of piezometric gradients, 
and they can be validated by the results given by models 
that have been implemented on a local scale. Biological 
markers that are or are not influenced by phreatic condi- 
tions and discharge fluctuations can also confirm these 
surface water/groundwater interactions. This novel ap- 
proach provides accurate maps of surface water/ground- 
water exchanges and demonstrates that hydraulics and 
biological metrics provide complementary information, 
despite the absence of significant correlation between 
hydrology and groundwater fauna, both in terms of 
spatial scale and the quality of information obtained. 
The results of this method will provide relevant hydro- 
logical and biological information on the complex interac- 
tions between groundwater and surface waters on the scale 
of the entire course of large rivers for decisions makers. 
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In terms of water resource management, the corres- 
pondence between hydraulics and biological indicators 
can help to determine sensitive areas where groundwater 
should be preserved. This is the case for the pumping 
well located in the southern part of the study area close 
to the Guiers River on the left bank of the floodplain. 
The maximum groundwater flow has been estimated to 
be 10.774 m7day over a 1-km distance on this left bank. 
Considering that the influence radius of the well is equal 
to 300 m, the extraction rate of the well (1.355 m /day) 
for supplying drinking water cannot be increased more 
than twice. In the case of groundwater overextraction, 
the exchange will be reversed, and the flow of the river 
can reach the pumping well. 

Future research will be required to validate the spatial 
correlation found between the metrics in another case 
study to take into account contrasting situations regarding 
latitude, climate, geological context, the biogeography of 
organisms and anthropogenic conditions. Similarly, tran- 
sient conditions remain to be studied because surface 
water/groundwater interactions can be modified during 
the fall low-flow period. Another field of research should 
be developed to predict the modification of these hy- 
draulic exchanges in the framework of climate changes. 
Indeed, the thermal modifications that are induced by glo- 
bal warming could have harsh repercussions on the water 
demand of human populations. For example, nuclear 
power plants would require increasing amounts of steno- 
thermal groundwater, and the groundwater demand for 
agriculture during summer periods would significantly 
increase. 

Additional file 



Additional file 1: Calculation of the trophic level of a wetland. The 

releve data (left table) is multiplied by the affinity of the species for each 
trophic class (center table) to obtain the right table, wherein the highest 
value (4.46) indicates the best estimation of the trophic level of the site. 
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